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Abstract 

V^ . A new approach to the solution of quasihnear nonehiptic first-order systems of 

CN I inhomogeneous PDEs in many dimensions is presented. It is based on a version of 

the conditional symmetry and Riemann invariant methods. We discuss in detail 

i-Q . the necessary and sufficient conditions for the existence of rank-2 and rank-3 

solutions expressible in terms of Riemann invariants. We perform the analysis 

'^ ' using the Cayley-Hamilton theorem for a certain algebraic system associated with 

Cd . the initial system. The problem of finding such solutions has been reduced to 

expanding a set of trace conditions on wave vectors and their profiles which are 
expressible in terms of Riemann invariants. A couple of theorems useful for the 
construction of such solutions are given. These theoretical considerations are 
^ , illustrated by the example of inhomogeneous equations of fluid dynamics which 

■^ I describe motion of an ideal fluid subjected to gravitational and Coriolis forces. 

OS ■ Several new rank-2 solutions are obtained. Some physical interpretation of these 

results is given. 

in 

o, 

^ : 1 Introduction 



^ 



In the last three decades, a number of useful extensions of the classical Lie approach to 
group-invariant solutions of PDEs have generated a great deal of interest and activity 
5^ ! in several fields of research, including mathematics, physics, chemistry and biology. 

In particular, the method of partially invariant solutions [131 IISI [HI [19], the non- 
classical method [21 [11], the conditional symmetry method [3 [6l [3 [T71 [18], the weak 
symmetry method [HIH] and the introduction of general "side conditions" or differential 
constraints [H] have been shown to play an essential role in several applications to 
nonlinear phenomena. In particular, the method of differential constraints incorporates 
all known methods for determining particular solutions of PDEs. They establish a direct 
connection between certain classes of solutions and a framework based on the theory 
of overdetermined nonlinear systems of PDEs. In an attempt to understand certain 
classes of solutions expressible in terms of Riemann invariants, we have introduced a 
specific version of the conditional symmetry method [6] . The basic idea was to analyze 
a Lie module of vector fields which are symmetries of an overdetermined system defined 
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by the initial system of equations and certain first-order differential constraints. It was 
shown [B] that this overdetermined system admits rank-fc solutions expressible in terms 
of Riemann invariants. The method of conditional symmetry has led to a large number 
of new results. 

This paper is a follow-up of our previous investigations which were performed in 
[21 El El E]- Here, we are using the Riemann invariants method for inhomogeneous 
nonelliptic quasilinear first-order systems of / PDEs in j9 = n + 1 independent variables 
of the form 

Ut + A\u)u^^ = B{u), (1.1) 

where x = {t,x^, . . . , x"') = {x^, x) G M"+^ and u = {u^, . . . , m*) G M"^ with the initial 
condition 

Xo = : u{0,x) = Uo{x) eR'^. 

As usual in mathematical physics, the space X of independent variables {xq, x) G 
W is called the physical space while the space of values of dependent variables u = 
{u^{x), . . . , u'^{x)) is denoted hy U (ZM.'^ and is called the hodograph space. The terms 
Uf, u^i denote the first-order partial derivatives of m with respect to t and x\ respectively. 
The matrices A^[u) are I x q matrix functions of u and B{u) is an /-component vector. 
Throughout this paper, we use the summation convention over repeated indices. All our 
considerations are local. It suffices to search for solutions defined on a neighborhood 
of X = 0. The solutions u = u{x) of the system (11. ip are identified with their graphs 
which are (n -|- l)-dimensional manifolds in the cartesian product X x U. 

In our previous analysis of ffist-order systems of the form (II. ip (see e.g. [B]), we 
considered homogeneous systems, i.e. systems with B{u) = 0, with coefficients de- 
pending only on the unknown functions u. The methodological approach proposed in 
these works was based on some generalization of the Riemann invariants method. A 
specific aspect of that approach is the presence of both algebraic and geometric points 
of view. A conversion of systems of PDEs into algebraic form was made possible by 
representing the so called integral elements as linear combinations of simple elements 
[6] associated with those vectors fields which generate characteristic curves in the X xU 
space. The introduction of those elements proved to be fruitful from the point of view 
of constructing the rank-k solutions in closed form by means of the Cayley-Hamilton 
theorem applied directly to the algebraic form of the system. Using a variant of the con- 
ditional symmetry method, we have shown that these solutions comprise among others 
multiple wave solutions as a superposition of two or more Riemann waves. The results 
obtained for the homogeneous systems were so promising that it seemed worthwhile 
to try to extend this approach and check its effectiveness for the case of nonelliptic 
inhomogeneous partial differential ffist-order systems. This is, in short, the objective 
of this paper. 

The paper is organized as follows. In Section 2, we investigate the group-invariance 
properties of rank-2 and rank-3 solutions for inhomogeneous systems (II. ip and obtain 
necessary and sufficient conditions for their existence. Section 3 illustrates the con- 
struction of such solutions for the inhomogeneous Euler equations describing a (3-1-1)- 



dimensional fluid flow under the influence of gravitational and Coriolis effects. Results 
and future perspectives are summarized in Section 4. 

2 Rank-/c solutions described by inhomogeneous sys- 
tems of PDEs 

A version of the conditional symmetry method [31 O [6] has been developed recently 
for homogeneous nonelliptic systems (11. ip . This method consists of supplementing the 
original system of PDEs f ll.ip with first-order differential constraints (DCs) for which 
a symmetry criterion for the given system of PDEs is identically satisfied. It turns out 
that under certain circumstances, the so-augmented system of PDEs admits a larger 
class of Lie symmetries than the original system of PDEs (11.11) . We now extend this 
method to inhomogeneous systems (II. ip and reformulate the task of constructing rank- 
k solutions expressible in terms of Riemann invariants in the language of the group 
analysis of differential equations. For this purpose, we require the solution u of PDEs 
(II. ip to be invariant under the family of commuting first-order differential operators 

d 

^a = C(w)^, a=l,...,p-k, k<p, (2.1) 

defined on X x U and satisfying the conditions 

CA| = 0, s = 0,...,fc-l, a = l,...,p-k. (2.2) 

Here, the inhomogeneous wave vector A° satisfies the rank condition 

rank (A°A\ B{u)) = rank (A°v4^) , (2.3) 

while the wave vectors A°, A^, . . . , A'^~^ are linearly independent vectors satisfying the 
dispersion relation for system (II. ip 

rank(A|A^) < /, s = l,...,k-l. 

The group-invariant solutions of the system (II. ip then consist of those functions u = 
f{x) which satisfy both the initial system (II. ip and a set of first-order differential 
constraints 

^a-g- = 0, a = l,...,p-k, a = l,...,q, (2.4) 

ensuring that the characteristics of the vector fields Xa are equal to zero. In this case, if 
the vector function u{x) is invariant under a set of {p — k) vector fields Xa for which the 
orthogonality property (12. 2p is satisfied, then the solution u{x) of (II. ip can be defined 
implicitly by the following set of relations between the variables m",x* and r** : 

M = /(rO,...,r'=-^), r'{x,u) = X'o{u)t + X'-{u)x\ s = 0, . . . , A; - 1. (2.5) 



Each function r^{x,u) is called the Riemann invariant associated with the vector A*. 

A vector field Xa defined on X x [/ is said to be a conditional symmetry of system 
(11 -ip if it is tangent to S = Sa H §q, where §a and E>q are submanifolds of the solution 
spaces defined by 

§A = {{x,u^^'^) ■.Ut + A\u)u^. =B{u)], 

§Q = { {x, «(!)) : ex. = 0, a = l,...,q, a = 1, . . . ,p - k] . 

An Abelian Lie algebra L spanned by the vector fields Xi, . . . , Xp^k is called a condi- 
tional symmetry algebra of the original system (11.11) if the following condition 



pI^^^Xa{ut + A\u)u,^-B{u)) 



0, a = 1, . . . ,p — k, 



is satisfied, where pr'^^^X^ is the first prolongation of X^. 

When studying solutions of type (12. Sp . it is convenient from the computational point 
of view to write system (II. ip in the form of a trace equation, 

TT[A''{u)du]=B^'{u), /i = l,...,/, (2.6) 

where A'^{u) are p x q matrix functions of u. The Jacobian matrix of relations (12. 5p 
can be expressed in matrix form either as 

du = «.) = f (x. - ivot + v.x^) ^) ' A, (2.7) 



or as 



where 



du =(!,-% ivot + V^x^)) '^A, (2.^ 



OJ (Of \ ^ -pqxfc \ — (\s\^Takxp „ I o\ 



* V*.;- • ' = (^)^^"- ". = (#j- ^' = "-■."■ 

and Xfc and Xg are the k x k and q x q identity matrices respectively. 

A solution of the form (12. 5 p is called a rank-fc solution if, in some open set of the 
origin x = 0, we can express u explicitly as a graph over the space of independent and 
dependent variables and the condition 

Ya.nk[du) = k 

holds. 

The condition of inversibility of the matrices appearing in equations (12. 7p and (12. 8p 

df df 



restricts the domain of existence of rank-fc solutions. However, since the inverse matrices 
are well defined at a; = 0, there exists a neighborhood of the origin in which det(Mj) 7^ 0, 
i = 1,2, allowing us to look for rank- A; solutions locally parametrized by equations (12. 5p . 
It should also be noted that because of the Weinstein-Aronzjain determinant relation 



det (Xfc - PQ) = det (X^ - QP), P E 



!)kxq 



Qe 



nqxk 



(2.9) 



which can be derived from the relation 



det 



^k 

Q 



p 



^k 

-Q 







det 



-Q 



Q 



p 



we have that det (Mi) = det (M2). Hence, the regions of validity for the expressions 
(12. 7p and (12. 8p are the same. Morever, (12. 9 p implies that det (Mi) is a polynomial of 
order s = min {k, q). Inserting (12. 7p into the original system (12. 6p . we obtain the result 



Tr 



A^ 



.u 



Xfc 






or, upon multiplication by det (Mi 



Tr 



df 



^^(u)adj ( Ik - {vot + Vix") -^ 



det 



Xt 



A =Bf'{u), 



{riot + riix') — 



B^'iu) = 0, 



(2.10) 

which is a polynomial expression of degree s in the independent variables (t, x'). For 
the system (I2.10p to be expressed only in terms of Riemann invariants, we must require 
that every coefficient of (t, x*) of this polynomial vanish identically. We then consider 
equation (I2.10p and all its partial derivatives up to order s with respect to the inde- 
pendent variables (t, x*) taken at the origin and set the resulting expressions to zero. 
The partial derivatives of the trace equation (I2.10p are obtained by making use of the 
Jacobi equations for the derivative of the determinant of any square matrix M = M{C,) 
depending on a parameter ^, 



|det[M] 
|adj[M] = 



Tr 



Tr 



adj[M] 



dM 



adj[M] 



dM 



d^ 



X-adj[M] 



dM 



(2.11) 



M~ 



A conditionally invariant solution of (II. ip under an Abelian algebra L spanned by vector 
fields Xi, . . . ,Xp_k of the form (12. ip is then obtained by solving the overdetermined 
system 



A': 



Tr 



^^(u)adj ( Xfc - {r]ot + n,x^) 
0, a = 1, . . . ,p — k. 



Therefore, in the case k 



dl 
dr 

a = 



A 



-det 
.,q. 



Ik - {vot + Vi^') 



dr 



2, we can devise the following proposition. 



Bf'(u) = 0, 



(2.12) 



Proposition 1. Consider a fixed set of linearly-independent wave vectors A° and X^ 
associated with a nondegenerate quasilinear nonelliptic first- order system of PDEs U.l\) 
in p independent variables and q dependent variables. This system admits a [p — 2)- 
dimensional conditional symmetry algebra L if and only if {p — 2) linearly independent 
vector fields 



Xa = Ca{u) 






A|C = 0, 



,p-2, 



s = 0,l, 



satisfy the conditions 
Tr 
Tr 
Tr 
det 



or 



^,9f^ df 



Or Or 



0, 



1,...,/, 



0, 



n, 



^8f f df df \ df ^ 

dr \ ^ dr ^ '^ dr ^ J dr 



0, 2 ^ J = 0, 



, n, 



df 



^ dr 



0, Vi 



9A| 
du° 



i)2xq 



(2.13a) 
(2.13b) 
(2.13c) 
(2.13d) 



and such that rank(9//(9r) = 2 on some neighborhood of a point {xo,uq) G S. The 
solution of ( li.ij) which is invariant under the Lie algebra L is precisely a rank-2 solution 
of the form l\2.5\) . 

Proof. It was shown in [6] that in the new coordinates MP x W^ 



x' 



r^{x, u),x'^ = r^(x, u),x^ = x^, ... 



^p 



X- 



x^u^ 



u 



m" 



u^ 



the vector fields Xa adopt the rectified form 



X - ^ 

O-f.a 



a = 3, . . . ,p. 



Hence, the symmetry criterion for G to be the symmetry group of the overdetermined 
system (12.121) requires that the vector fields Xa of G satisfy 

XaiA')=0, 



whenever equations (12.121) hold. Thus the symmetry criterion applied to the invariance 
conditions (12. 4p is identically equal to zero. After applying the symmetry criterion to 
the system (I2.10p in new coordinates and taking into account the conditions (I2.13p . we 
obtain the equations which are identically satisfied. 

The converse is also true. The requirement that the system (11.11) be nondegenerate 
means that it is locally solvable and is of maximal rank at every point {xo,Uq 



(l)^ 



§. Therefore 



the infinitesimal symmetry criterion is a necessary and sufficient 



condition for the existence of the symmetry group G of the overdetermined system 
(12 ■121) . Using relations f l2.1ip . we show that equations f l2.13p must hold. Equations 
fl2.13ap and fl2.13bp are easily obtained by considering the system fl2.10p and all its 
first-order partial derivatives with respect to t and x*. The second-order derivatives 
with respect to x* and x^ provide us with the condition 



Tr 



\Bf f df df \ df ^ 

dr \ ' dr '^ "^ dr ^ j dr 



0, 2,J=0, 



which reduces when i = j to 

Tr 



dr dr dr 



0, i = 0, 



n. 



(2.14) 



Considering equations (I2.13bp . we can write (I2.14p as 



Tr 



dr dr dr 



Tr 
Tr 



'^,dl ^dl dl^ 
dr dr dr 



Tr 



' .u9f df f df ^ 
dr ^ dr \ ^ dr 



dr 

\ — 
dr 



Tr 
X2^ A 



Audf df ^ 
dr dr 



(2.15) 



According to the Cayley-Hamilton theorem, each matrix [rji-^^ ^^^"^ satifies 



dl 
' dr 



Tr 



' dr 



dr 



det 



3/ 
' dr 



Xo 



allowing us to write (I2.15P as 



Tr 



■ df df df^ 
dr dr dr 



det 



dr 



Tr 



dr 



Bf^det 



■ df 
dr 



by (I2.13ap . As a consequence of the assumption B'^ ^ 0, we must have that 



det 



^ dr 



0. 



That ends the proof, since the solutions of the overdetermined system (I2.12p are invari- 
ant under the algebra L generated by the (p — k) vector fields Xi, . . . , Xp_2. The invari- 
ants of the group G of such vector fields are provided by the functions {r°, r^ , u^ , . . . , u'^} . 
So the general rank-A; solution of (II. ip takes the form (12. 5p . 



D 



Note that in the case q = k = 2, equations (]2.13dp imply that 



det 



dr 



det [rji] det 



df 
dr 



0. 



When det \-^'\ = 0, the considered solution reduces to a rank-1 solution. Consequently, 
when q = 2, rank-2 solutions exist only when det [rji] = for i = 0, . . . ,n. Note that 
in every case, equations (I2.13dp can be satisfied, for example, by choosing one of the 
vectors A* to be constant. 

Analogously, we obtain the following conditions in the case k = 3. 

Proposition 2. Consider a fixed set of linearly-independent wave vectors A°, A^ and A^ 
associated with a nondegenerate quasilinear nonelliptic first- order system of PDEs U.l\) 
in p independent variables and q dependent variables. This system admits a {p — 3)- 
dimensional conditional symmetry algebra L if and only if {p — 3) linearly independent 
vector fields 



Xa = Cain) 



_d_ 



satisfy the conditions 



Tr 
Tr 
Tr 
Tr 



B^, 



or 

'a'^%. ^a 

dr ^ dr 
' df f df 



/^ 



m: = 0, 



J, 



,p 



0,1,2, 



^'^ > ^- 



.df 



0, ii = 0, . . . ,n 
df \df 



"1^2 + Vi2-Q^Vh 



dr 

df df 



dr 



A 



0, ii,i2 = 0, 



(2.16a) 
(2.16b) 
(2.16c) 



or V or or or or or or 



df df _ „, , _ , 



dl di 

dr *^ dr 



df df \ df 



det 



' dr 



0, m 



*3 Q^ ■"! Qj, 

d\i 



du"^ 



3xq 



(2.16d) 
(2.16e) 



and such that rank(9//9r) = 3 on some neighborhood of a point (xq, uq) G S, where the 
indices ii, i2, is take the values 0, . . . , n and are not all equal. 

Proof. The proof for sufficiency follows the same steps as in the case k = 2. We now 
show the necessity of equations (12.161) . As previously, equations (I2.16b -d) represent the 
first, second and third-order partial derivatives of expression (I2.10p taken at the origin, 
respectively. When ii = i2 = is = i, equations (I2.16dl) reduce to 



Tr 



df df df 9/ ■ 
dr dr dr dr 



0. 



For any 3 by 3 matrix M, the Leverrier-Faddeev algorithm and the Cayley-Hamilton 
theorem allow us to write 

M^ - Tr [M]M^ + - (Tr [Mf - Tr [M^] ) M - det[M]X3 = 0. (2.17) 



Tr 



Hence, taking into account equations fl2.16p and (12.171) . we get 



dr dr dr dr 



Tr 



Vi 



dr 



Tr 



dr ^ dr * dr 



dr \ dr J \ dr 



^ dr 



X.]\ 




r]i—- + det 
dr 



dr dr 



+ det 



■ df- 
dr 

m— 



Tr 



A 



dr 



which is equal to zero if and only if the relations 



det 



Vi 



dl 
dr 



0, 



0, 



n 



hold. Again, the solutions of the overdetermined system (12.121) are invariant under the 
algebra of vector fields Xi, . . . , Xp_^ defined by (12.11) . (12. 2p and (I2.16p . The invariants 
of the group G of such vector fields are provided by the functions {r°, r^ ^r"^, u^, . . . , u'^}, 
therefore the rank-fc solution of (11.11) is of the form (12. 5p . D 



3 Applications in fluid dynamics 

Now we present some examples which illustrate the theoretical considerations presented 
in Section 2. We consider classical equations of hydrodynamics describing a motion in 
a fluid medium when the gravitational force g and Coriolis force Q x v occur. We 
restrict ourselves to the equations of a one-component nonviscous fluid flow. Under 
these assumptions, our equations are of the type (II. ip . The matrix form of these 
equations of hydrodynamics in a noninertial coordinates system is 



A°{u)ut + A^{u)u,^ + A'^{u)uy + A^{u)u^ = B{u) 



(3.1) 



where 



A'^iu) 



Bin) 



/ 1 \ 
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1 

1 

\ -Kp/p 1 / 

f gi- (^2^^3 - ^3^2) \ 



V / 



AHu) 



u 



• P,P) 



Vi 
V, 







Vi 








Vi 




S^l/P \ 
5i2/p 

Sis/p 









— KpVi 


Ip 


V^ J 


GM^ 


Sij 


={: 


i 
i 


= 3 



1,2,3, 



Here we treat the physical space X C M^ as classical space-time, each of its points having 
coordinates (t, x) = {t, x, y, z) and the space of unknown functions (i.e. the hodograph 
space) f/ C M^ having coordinates {v,p,p). The vector ^2 denotes the angular velocity 
while K is related to the adiabatic exponent 7 of the fluid by the relation k = 2(7 — 1)~^. 
For the homogeneous system (11.11) . each wave vector A can be written in the form 
A = (Ao,A), where A = (Ai,A2,A3) denotes a direction of wave propagation and the 
eigenvalue Aq of the matrix (AjA*) is a phase velocity of the considered wave. The 
dispersion relation for the homogeneous hydrodynamic system of equations (13. ip takes 
the form 



(Ao + ^-A)= 



(Ao + A 



— lAp 
P 



0. 



(3.2) 



Solving equation (13.21) . we obtain two types of wave vectors, namely 



|A|-t;-A,A , £ = ±1. 



i) Entropic : A = [—X-v,X), 

\ 1/2 
ii) Acoustic : A^" = i e { — 



Algebraic equations deflning inhomogeneous wave vectors which satisfy equation (12. 3p 
are of the form 



i) Entropic : Xeo 



ii) Acoustic : A^o = ( ^ ( — 



iii) Hydrodynamic : A/^o 



-V ■ g,g — Q X V 

1/2 



|A| — -u • A, A , , 



Aq, a 



Ao + w- A ^0, 



±1 



Kp 



P 



1/2 



lAI 



Several classes of wave solutions of the hydrodynamic system (13. ip have been obtained 
via the generahzed method of characteristics [U [10]. Applying the CSM to this sys- 
tem allows us to compare the effectiveness of these two approaches. Here we restrict 
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ourselves to the consideration of superpositions of waves which are admissible by the 
inhomogeneous system fl3.ip . We are interested only in rank-2 solutions which can be 
written in Riemann invariants. The results of such superpositions are set in Table [H 
where the signs (+) and (-) indicate when superposition does and does not occur, re- 
spectively. We will now show for illustration some physically interesting solutions from 
among the ones obtained. 





E 


Ae 


E^ 


+ 


+ 


A^s 


+ 


— 


H' 


+ 


+ 



Table 1: The existence (+) or absence (-) of nonlinear superposition of waves admitted 
by the system (13.11) . 



Different possibilities of existence of solutions in Riemann invariants for different 
combinations of wave vectors (homogeneous and inhomogeneous) are denoted by sub- 
systems E^E^ E^A^, A^^E, A^Ag, H^E, H^A^. Moreover, we denote by r° a Riemann 
invariant associated with wave vectors A^, A^^ , A^ and by r^ a Riemann invariant as- 
sociated with A^, A^^ 

Now we present some examples which illustrate the theoretical considerations pre- 
sented in Section 2. 

Subsystem {E^, E} : We look for solutions of the form (12. 5p defined by 



v{r ,r ), p = p{r ,r ), p = p{r ,r ), 

■ -{v- g)t+ {gi + Q3V2 - 112^^3)3; + {92 + ^iV3 - ^3Vi)y + {93 + ^2t^i - ^1^^2)2;, 
X'-v)t + Xlx + \\y + A^;., A^ = (Ai, A^, A^), 



(3.3) 



where the A] are allowed to depend on the unknown functions. To satisfy equations 
(I2.13d[) . we require that the coefficients of t,x,y,z in r° be constant, say r° = cot + 
C\X + c^y + c-^z. To allow for a rank-2 solution, the following relations must then hold 



^■fi = 0. 



V\ 



g'2 - C2 + VL\V3 



V2 



ci- gi + ^2f 3 

^3 



implying that this type of solution exists only when the rotation axis of the system 
is perpendicular to the constant gravitational force. Solving (12.131) . we obtain the 
following solution 



Vl 



Co^2 



91 



Qo 



"3 9i 9i93''h ^h 93^h ^h ^h 



(3.4) 



p = PrO, P = p{r^), 
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where p{r ), f3(r , r ) are arbitrary functions of the Riemann invariants 
r° = - {gst + Q2X - Qiy) , g-n = 0, 



93 



>^l9i - >^\92 , CqXI 



n. 



93 



+ ^^^ ] t + X\x + Xly + Xlz, 



(3.5) 



and A^ is a constant vector satisfying A^ • i7 = 0. The solution defined by ( 13 ■4p and 
(13. 5p represents a double traveling wave with constant velocities. 

Subsystem {E^, A^} : We look for solutions of the form (12. 5p defined by 

t+{gi+ ^3172 - ^2^^3)3; + {92 + ^iV-i - ^3Vl)y + {93 + ^2^^1 - ^lV2)z 



r = — [V 



np 



1/2 



where the A] are allowed to depend on the unknown functions. The solution in this 
case exists only when k, = 1. Following the same process as with subsystem {Ei,E}, 
the rank-2 solution exists when g ■ Q = and is given by 



Vl = C^ 



\n\ 



^iB(r'] 



V2 



QiQs 



'' = n; 



cifii - giili + ^2 \ C2 - 92 + ^3 \ Ci 



M 



Q,B{r'] 



C2- g2 + ^3 Ci 



M 



Q.Bir'] 



Here B{r^) is an arbitrary function, the Riemann invariants are given by 

C2gi-Cig2, , , 01^1 + 02^2 
r = 1 + C\X + C2y ;: 2, 



^3 

A 



n. 



\n\ 



{X^ ■ n)B{r^) - xlCi + Va\x^\ 



, Xln^jgi - ci) - (A^l^a + Xlns){c2 - g2 + €1^3) ^ , ,1 ,,1 ,,1 

H — — j I + A^X + A2y + ^3^? 



(3.6) 



(3.7) 



and A^ is colinear with fi, A^ x fi = 0. Note that this solution admits the gradient 
catastrophe for a certain time Tq that we estimate by considering a first-order taylor 
expansion for B[r^) ~ i?o + Bir^. We obtain that the first-order partial derivative of 
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(13 .yp with respect to t goes to infinity at time 



Tn 



BiVA(x^-n 



(3i 



A physically interesting solution can be built by selecting B{r^) as the Jacobi elliptic 
function sn(r^,m), where m G [0,1] is the modulus. For this choice, solution f l3.6p . 
becomes 



vi=Ci 



V2 



\n\ 



Qi sn(r^, k), 



^1^3 



^3 



fil 



ciQi - giQi + ^2 \ C2 - 92 + ^3 { Ci 
A 



\n\ 



fii sn(r^,m) 



(3.9) 



C2- 92 + ^3 Cl 



\n\ 



rii sn(r^,m) 



while the invariants adopt the implicit form 

C29l-Ci92^ , , CiQi+C2^2 
r = ;: 1 + C\X + C2y — Z, 



A 



n. 



\ft\ 



(A^ • n)sn{r\ m) - X\Ci + VA\\^\ 



(3.10) 



, A^^i(^i - Cl) - (A^^2 + A^1]3)(C2 -92 + Cias) u , M ^xi ^xi 
H — — ) t + A^x + A2I/ + A^z. 



Figure [T] illustrates the evolution of the density function p{r^,r^) in (13. 6 p near the 
gradient catastrophe when the function B{r^) is chosen to be one the Jacobi elliptic 
functions sn(r^, 1/2), cn(r^, 1/2), dn(r^, 1/2) for a certain choice of the parameters. It 
should be noted that for solution (13. 9p . the gradient catastrophe can be estimated to 
take place, according to formula (13. 8p . around the time 



To 



101 



y/A{X^ ■ Q) 
Subsystem {A^, E} : We now look for solutions of the form (12. 5p defined by 



-\\^-v + e{'^ 



1/2 



|A°| t + A?s + A°y + A°2, 
r' = -(XK^y + Xlx + Xly + A^^, A^ = (Xl A^, A^), 



5 = 0,1, 
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2 1 6 



-1 _2 1 



-1 -2 




2 1 6 



-1 _2 1 



-1 -2 






2 1 



-1 _2 1 



-1 -2 




2 1 



-1 _2 1 



-1 -2 



Figure 1: Density distribution for the solution (13. 6p at time t = and near the gradient 
catastrophe for the eUiptic functions B{r^) = sn(r^, i), cn(r^, |), dn(r^, |) . 



where the Aj are allowed to depend on the unknown functions. In this case, the solution 
exists only when A'' = ei|A°|r2, \Q\ = 1 and g ■ Q = 0. Under these circumstances, the 
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rank-2 solution takes the form 



P = Po, P = Po, Po,Po^ 
1 

^3 



n-iEi 



^VPof^/p - ^l(^l^^l + ^2t^2) 



Co 



where the functions Vi, V2 satisfy the equations 

f 1 \ I di d2 \ I vi 



with the constants being given by 



V2 I \ df 



(3.11) 



di 



^1^2 



eQs^/poK/po 
ni^ + Qs^ 



, d2 



^2' + ^3' 



(^5 



eVL^^JpQK\/po' 

VL1VL2 

^^sVpokJpo' 



d?. 



£^^3a/?WPo 
eVL^^Jpqk/pq 



The eigenvalues p± and eigenvectors V± of matrix ( ^ ^ ) ^'^'e given by 



Po 



p± = ±eij , V± 

PqK 



-1 



The solution of equations (13. lip is then given by 

vi = Fi(r^) cos \e^/po/Kpor°) + F2{r^) sin \e^/po/Kpor°) + ci, 

Fi{r^) cos (^e^/po/npor'^] + F2{r^) sin \e^/po/Kpor^^ + cij 



V2 



+ 



+ 



gl^3 
£(fi2 + ^2 

£^2y/poK/po - C0Q2 - eifi'1^2^3 



F2(r^) cos ( £:Va)/Vo'^° ) " ^i(^^) sin (e^/po/Kpt 



Jor 



ei{ni + ni) 



where Cq G M, Fi(r^), F2(r^) are arbitrary functions of the Riemann invariant r^ and 
Ci G M is given by 



Cl 



£i(g'2^i + ^^(5^2 + g'l^i)) - C0Q1Q3 + eQiQs^ypon/po 



The Riemann invariants then adopt the explicit form 

r° = Cot + El {Qix + VL2y + ^3-2;) , 

r^ = £1 [e^JpQn/pQ -cojt- Qix - ^2!/ - ^3^- 
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This solution models the superposition of two traveling waves with the same direction 
but different phase velocities. 

Subsystem {H^E} : We finally consider the case of the propagation of an entropic 
wave on an hydrodynamic state. Since the required computations are very involved, we 
first assume, without loss of generality, that the coordinate system is chosen in such a 
way that the vector fl is given by (0, 0, 1). Using equations fl2.13p . one can show that a 
laminar solution of equations (13.11) exists when the wave vectors A*^ and A^ are colinear 
and is given by 



Vl 



5-2 H exp 

C2 



-W 



V2 = exp 



-W 



C2 (co + cig2 - C2gi) 



^1, 



C2{co + cig2 -C2gi), 
V3= I ^ — dr° + F2(ri), 

Co + CiVi + C2V2 

C2 [Co + Cig2 - 025-1) 
where the Riemann invariants satisfy the relations 



(3.12) 



Cot + CiX + C2y, C2 7^ 0, 



2 I 2 

Cl +C2 



1, ^ = (0,0,1) 



cig2 - C2gi + exp 



-w 



t + cix + C2y. 



(3.13) 



^2 (co + cig2 - C2gi] 

Here, the functions Fi(r^), F2{r^) are arbitrary functions and p{r^) is an arbitrary pos- 
itive function of the Riemann invariant r^. The function W{(f)) stands for the principal 
branch of the Lambert W function which is defined by the relation 



W (0) exp [W (0)] = 4>. 



(3.14) 



The Lambert function W{(j)) is positive and increasing when its argument is greater 
than and negative only for — e~^ < < 0. Moreover, it is unbounded when tends 
to +00. Hence, the velocity components fi,f2 are bounded for positive values of 



C2(co + cig2 - C2gi) 

Also, the solution is perturbed in accordance with the nature of the arbitrary function 
Fi(r^). This situation is illustrated in figure |2] for the case of a periodic perturbation. 



smr 



4 Final remarks 

The main result of this paper consists of extending the Riemann invariant method to 
inhomogeneous quasilinear systems with coefficients depending only on dependent vari- 
ables u. The conditional symmetry method, together with the Cayley-Hamilton idea. 
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Figure 2: Velocity vector field V2 for solution ( 13.12^ at times t = and t = 1 when 
Co > and Fi(r^) = sinr^. 

was adapted to partial differential equations in such a way as to allow the applicability 
of the Riemann invariants method. The conditions fl2.13p and fl2.16p are necessary and 
sufficient for the existence of rank 2 and 3 solutions of the system (11 .ip , respectively. In 
the proof of theorems 1 and 2, we obtain that the gradient of the matrix functions M^ 
or M2^ becomes infinite for certain values of x° (the solution and its derivatives do 
not remain bounded). This means that we are dealing with a gradient catastrophe and 
some discontinuities can arise such as, for example, shock waves. Thus, the existence 
of Riemann invariants implies that it is not possible to introduce the concept of a weak 
solution in the broad sense. 

These theoretical considerations have been illustrated by the examples of inhomoge- 
neous equations of fluid dynamics in the presence of gravitational and Coriolis external 
forces. We were able to construct new classes of rank-2 solutions expressible in terms of 
Riemann invariants in the cases E^E, E^A^, A'^E, H^E, H'^A^. Some of these solutions 
admits a gradient catastrophe. This results from the fact that the first derivatives of 
Riemann invariants become infinite after some finite time. These solutions, in their 
general form, possess some degree of freedom, that is, depend on one or two arbitrary 
functions of one variable (Riemann invariant) depending on the case. This arbitrariness 
allows one to change the geometrical properties of the governed fluid flow in such a way 
as to displace the singularities to a sufficient extent. This fact is of some significance 
since even for arbitrary smooth and sufficiently small initial data at t = to^ the mag- 
nitude of the first derivatives of the Riemann invariants become unbounded in some 
finite time T. This time can be estimated for each solution. 
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